Dynamics of a self gravitating light-like matter shell with spherical symmetry 
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Complete hamiltonian description of the physical system composed of a null matter shell inter- 
acting with the gravitational field is provided. In spherically symmetric case, the phase space of 
the system is effectively reduced with respect to the Gauss-Codazzi constraints. The Hamiltonian 
of the system (numerically equal to the value of the ADM mass) is explicitly calculated in terms of 
the „true degrees of freedom", i.e. as a function on the reduced phase space. Geometric interpreta- 
tion of the momentum canonically conjugate to the shell's radius is given. Transformation between 
different time parameterizations of the shell is discussed. 
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I. INTRODUCTION 

Null matter shell may be considered as a singularity 
of the space-time geometry consisting of two different, 
smooth 4-geometries tailored together along a hypersur- 
face S, where "tailoring" means that the (properly un- 
destood) induced geometry of S is continuous (cf. 
Singularity arises because the four-dimensional connec- 
tion coefficients F^^ may be discontinuous on S. Hence, 
the singular part of the Einstein curvature tensor den- 
sity of such a spacetime contains derivatives of those dis- 
continuities and may be nicely described in the sense of 
distributions as Q\ = G°(,<5s, where ds is a Dirac delta 
distribution concentrated on S. In paper Q dynamics of 
the composed "null-matter-shell -I- gravity" system was 
derived from first principles, in both Hamiltonian and 
Lagrangian approach, and for a general null matter. 

For this purpose an appropriate Lagrangian L — 



grav 



'-^matter was uscd. For its gravitational part igrav, 
R was taken, with 4-dimensional 



where i?ioe; is the 



L 

the quantity j^y/\g\ 
curvature scalar R 
regular part of the scalar curvature and \/]g\Rsing = 
— sing(CJ) = — G^'^5/ii/<^s is its singular part. (For cal- 
culational purposes, an adapted coordinate system may 
always be chosen such that S := {x^ = const.} and, con- 
sequently, we have 6s = S{x^).) 

In case of a spherically symmetric shell this formula- 
tion leads to a simple Hamiltonian system. Complete 
hamiltonian description of a massive dust shell was pro- 
posed in Q and Here, we are going to derive similar 
results for the case of a null shell. Similarly as was done 
in [i*], spherical symmetry condition is imposed on the 
Hamiltonian level only. (On the other hand, one could 
also impose this symmetry already on the level of varia- 
tional (Lagrangian) formulation, as was done in Q.) As 
a result, we obtain the Hamiltonian dynamics of the sys- 
tem described in terms of gauge-independent variables 
and their conjugate momenta. 

Fundamental difference between the null and the mas- 
sive case consists in the fact, that there are no massive 
shells without matter. On the other hand, gravitational 



"shock waves", i.e. null shells without any matter, are 
perfectly allowed by the theory. Nevertheless, our theory 
admits also null-like matter, which couples consistently 
to gravity. We have analyzed several models of a La- 
grangian function describing null matter. In most cases 
they couple to an isolated horizon only, which is rather 
a trivial case. However, interesting and self-consistent 
models have been found (see Appendix E]) . which lead to 
non trivial dynamical systems. 

The space of initial data of a "matter -I- gravity" sys- 
tem can be parameterized by the following space of func- 
tions: 

7':={(g,,,P",z^,PK)} . 

To remind, gki is a three-dimensional metrics on a space- 
like Cauchy surface, and P'^' are appropriate ADM (Q) 
momenta describing external curvature of this surface. 
Moreover, describe configuration variables of the mat- 
ter and pk, their conjugate momenta. We limit ourselves 
to the topologically trivial case C ~ R"^ and assume that 
the geometry of C is asymptotically flat in infinity. 

The following phase space is equipped with the canon- 
ical pre-symplectic structure: 



1 

16^ 



+ J {SpK A Sz'^) (fx . 



(1.1) 



Assuming that the matter is concentrated on the shell S 
exclusively, we obtain: 



1 

16^ 



{5P''\x)h5gki{x)]d' 



+ / {SpK A Sz'^) (fx 
Jens 



(1.2) 



The structure of the paper is following. In Section ITTl 
we analyze the above pre-symplectic form and reduce it 
with respect to both the spherical symmetry and the con- 
straint equations fulfilled by ADM data. It turns out that 
the radial component of constraint equations gives us an 
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"equation of state" , which has been previously postulated 
in paper Q. Exact solution of Gauss-Codazzi constraint 
equations are obtained in Section lTTTl Using these results, 
both the gravitational part (Section IIV Ap and the mat- 
ter part (Section llV Cp of the symplectic form is reduced. 
Then, an exact form of the Hamiltonian of the system is 
derived, which encodes its equations of motion. These 
equations are solved with respect of the Schwarzschild 
time in Section |Vl Then we perform transition from 
Schwarzschild to Minkowski time (Section IVip . and solve 
resulting equations of motion. In Appendix [Xl we discuss 
different models for the null matter. 

This paper is a continuation of a serie of papers ^ and 
Most of the results described here have been included 
in the Ph.D. thesis ^ of one of us (E.G.). 



II. SPHERICAL SYMMETRY 

Spherical symmetry of the system implies the existence 
of spherical coordinates {x^,x'^,x^) = {9,(p,r) on the 
Gauchy surface C, such that initial data {gki, P'''' , ,pk) 
are invariant with respect to rotations. Moreover, we 
may choose the radial coordinate in such a way that it 
is constant on the shell (i.e the history of the shell is 
given by equation r = = — const). This is not a 
physical assumption, but merely a choice of a gauge, as 
legal as any other gauge fixing condition (e.g. conditions 
used in papers [1, [l^, [13 ) • It simplifies considerably our 
calculations. 

Another gauge condition, which we use in this paper 
is the continuity of all the 10 components of the met- 
ric. In many approaches only the geometry of the shell 
is supposed to be continuous across the shell. Our condi- 
tion is stronger. Again, this is not a physical restriction 
imposed on the theory but merely a gauge fixing con- 
dition which allows us to use theory of distributions to 
calculate the (singular!) curvature tensor. In particular, 
our condition implies equivalence of the afFine structure 
of null geodesies on the shell, as seen from the external 
and the internal geometry. Without our gauge condition, 
this equivalence must be imposed separately (see [11 for 
further discussion of the role of this condition). 

Due to spherical symmetry, the two parts of the ge- 
ometry, which we taylor together through the shell, must 
be Schwarzschild. Moreover, due to topological trivial- 
ity of the Gauchy surface C, the internal part must be 
Minkowski. 

We assume that asymptotically, for r — > oo and r —f 0, 
the radial coordinate is "asymptotically flat" , i.e. equals 
to the Schwarzschild radius outside and the Minkowski 
radius inside the shell. 

The three-dimensional metrics on C must, therefore. 



be of the following form: 





" 1(1^) lAB 


" 


gki 











n^(r) _ 



where I and n are functions of radial coordinate r, and 
lAB (where A,B — 1,2 label angular coordinates) de- 
notes standard metrics on the unit sphere: 

lAB^{l .\\ . (2.2) 
\ sm j 

Spherical symmetry implies also that the components P^^ 
of ADM momenta assume the following form 



n[r) 

= , (2.3) 

where u and / are again functions of the radial coordinate 
r. These functions are piecewise smooth outside of the 
surface of the sphere given by equation r — C^^ whereas 
the metric coefhcients n and I are supposed to be also 
continuous at r = ^. 

Einstein equations imply constraints which must be 
satisfied by the above data, namely the Gauss-Godazzi 
equations for the components 5*^^ of the Einstein tensor 
density. Standard decomposition of Cy''^ into the spatial 
(tangent to Vt) part and the time- like (normal to Yt) part 
gives us respectively: 

Q\ = -Pi\k , (2.4) 

1Q\n^ = -sg R +s (^P^'P^i - ^P^^ ^ , (2.5) 

where n is the future orthonormal vector to Gauchy sur- 
face Vt and s is a sign of g'^". In paper [13] (see also [1|) 
we have shown how to decompose these equations into 
the regular and the singular (proportional to Dirac delta 
distribution on St) parts. The regular part of the vector 
constraint reads: 

reg(P/V)=0' (2-6) 

whereas the regular part of the scalar constraint reduces 
to: 

reg (^g R - {p'^'Pti - \p^^ ^) = ° ' ^^'^^ 

The singular part of divergence of the ADM momentum 
Pki contains derivatives in the direction of x^ : 

sing(Pi'=|fc) = sing(a3Pi') = S{x^)[Pi^] , 
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so the right-hand-side of (|2.4p has the form 

Pi''\k = [Pi']Six^) ■ (2.8) 

where the bracket denotes the jump of Pi^ between the 
two sides of the singular surface. In paper we have 
proved that components P*^' are regular, whereas the 
singular part of the three-dimensional curvature scalar 
reduces to the jump of the external curvature k of the 
two-dimensional surface St C Vt across the shell: 

(3), 



sing(5 R) = 2g^[k]S{x'') = 2^l[k]Six^) 
Hence, the right-hand-side of l|2.5p has the form: 



sg R -,s (^P'^'Pki - ip^^ i = 2^l[k]dix') . (2.9) 

Due to Einstein equations, these objects must match the 
matter energy-momentum tensor. A priori, there are se- 
rious difficulties concerning definition of such an object 
in case of a matter living exclusively on the null surface 
S. We have shown in paper Q that such a quantity 
may be consistently defined. It is a three-dimensional 
tensor-density denoted by r°fc (one three-dimensional in- 
dex a = 0, 1, 2 up and one index down!). We stress that, 
due to the degeneracy (-^det gat = 0) of the metric tensor 
gab on 5, it is impossible to find any tensor representa- 
tion of this quantity. Also the, otherwise trivial, "rising 
of indices" is forbidden, whereas the "lowering" becomes 
a non-invertible (loosing information) procedure! 

Naively, one could expect Gauss-Coddazzi constraints 
in the form: 



We stress, however, that there is no way to define the 
right-hand side as components of any well defined four- 
dimensional object! In our notation, where x^ is constant 
on 5, only t'^i, makes sense! Fortunately, we have shown 
in 0] (see also [l2|) that, due to the null character of S, 
the singular part of constraint equations contains only 
three conditions. Indeed, the orthogonal (to S) part of 
the constraint equation coincides with one of the tangent 
parts of these equations, namely the component along 
the null vector on S. (The fourth constraint, existing in 
a non-degenerate case, is replaced here by the degeneracy 
condition det gab for the metric on S.) The tangent (to S) 
part of G^b splits into the two-dimensional part tangent 
to St and the transversal part (along null rays). The first 



one gives us: 



[P'b] = -8nr% . 



(2.10) 



In spherically symmetric case both left- and right-hand- 
side must vanish and the above constraint is fulfilled au- 
tomatically. The remaining null tangent part of Einstein 
equations reduces, as we have shown in [3|, to the follow- 
ing constraint: 




i/det gAsk 



. 



(2.11) 



In our notation it can be rewritten in the form, which 
has been postulated in paper Q: 



- s- + e = 
n 

where we have denoted: 
-^,:=i^sing (P33|^)=±[P33]^ 



(2.12) 



(2.13) 



e:^^sing R -(P^^ - Ip^)^) ^ ^^VVm 



(2.14) 



The "equation of state" (|2.12p is a consequence of "nuU- 
ness" of matter and does not come from its internal prop- 
erties. 

Equations (prUl) . (PTTjl . together with ([M]) and 
provide the complete set of constraints fulfilled by the 
initial data {l{r),n{r),u{r), f{r)) of the theory. 

To reduce the phase space with respect to these con- 
straint we will proceed as in a paper [j], but with some 
modification. Equation (|2.13p . together with the vec- 
tor constraint (|2.6p . can be written in terms of momenta 



O-f-^-l' , (2.15) 
2 n 



ins— A(r 
n 



where prime ,,/" denotes the radial derivative d/dr and 
A is the one-dimensional Dirac delta. Equation (j2.14p . 
together with (|2.7p . may be rewritten in an analogous 
way: 



-87re(p,;,n) A(r-C) 



-n--— — f^~-fu . 

A In Al 2-' 

(2.16) 

The above constraints generate two dimensional group 
of space-time reparameterizations, where variables {t, r) 
may be replaced by any other arbitrary variables (t,r), 
preserving spherical symmetry of our system. Gauge 
transformations enter as degeneracy directions of sym- 
plectic structure Vsym, obtained by restriction of the 
form Q from V to the space of spherically symmetric 
data which we denote by Vsym- In order to calculate this 
restriction let us consider equation (|1.2p in case our spe- 
cific case and integrate over angular momenta. Hence we 
obtain the following symmetric structure in Vsym- 



n = I [ ]su{r) A dl{r) + ^^/(r) A Sn{r) ] dr 



A7t6{PkI{(:)) a Sz^ 



(2.17) 



where we denote: pK ='- \/det gAsPx- Moreover, for 
purely calculational reasons it is convenient to write 2- 
form Q as an exterior derivative f2 = SQ of the following 
1-form: 



e 



1 



u{r)dl{r) n{r)5J{r) dr 



■ATTliQPKSz 



K 



(2.18) 
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As already mentioned, Cauchy surface C is tailored 
from a piece a 3-surface in the flat Minkowski spacetime 
for r < ^ and a piece of a 3-surface in the Schwarzschild 
space-time for r > 0. 



III. SOLUTION OF CONSTRAINT EQUATIONS 

In order to solve constraint equation we have to impose 
gauge condition which enables us to fix uniquely the time 
coordinate. For purely technical reasons we start with 
the family of /?-gauge conditions, used for purposes of 
the proof of positivity of the ADM mass in paper p^ . 
The condition reads: l3P^^g33 + P^^ gAB = 0, where (3 is 
a fixed constant. This means, in terms of our variables: 



u 
n 



~/3 



Assume that /3 < — 1. Inserting the above relation to the 
vector constraint (|2.15p we see that, outside the shell, our 
initial data must fulfill the following equation 



(3,2) 



This implies that the function log(/ l^) is constant ev- 
erywhere outside of the shell. Hence we have: 




dla r > C 
dla r < C 



(3.3) 



The jump — is determined by the constraint 
equation l|2.15p - i.e. the only singular term on the right 
hand side of that equation is produced by the jump of 



/ and is equal to {A+ - A_)(Z(C))" 
physical radius of the shell by 



R 



2 A(r — C). Denote 



and normalized momentum radial component by U 
Equation (|2.15p implies: 



(3.4) 



(3.5) 



To solve completely constraint equation, the boundary 
conditions must be imposed, i.e. for r ^ oo initial data 
must be asymptotically flat, and for r — > must be reg- 
ular. This implies that / must behave like both at 
zero and at inflnity, whereas / must vanish at inflnity. 
Consequently, for /3 < — 1, we obtain A+ = 0: 



/ = 







for r > ( 



2UR^+H-'2 for r < C 



(3.6) 



Now we are going solve the scalar constraint equation 
(|2.16p . Using the gauge condition (|3.ip and the vector 
constraint (|3.3|) . we obtain the following equation, which 
is valid outside and inside of the shell: 







4 In 



1 + 2/? ,2 



Consider the following quantity: 

V 



k := 



We see that (|3.7I) may be rewritten as follows 



^k'-l+'-^Al 



(3 1) equivalently, 



4/1+/3 



Al 



(3.7) 



(3.8) 



(3.9) 



(3.10) 



This means that the function in brackets must be piece- 
wise constant. Regularity of the metric at r = implies 
that for small r the function I shall behave like n^r^. This 
implies that k'^ behaves like 4-\/I and, consequently, the 
corresponding internal constant vanishes. Denoting the 
remaining external constant by ,,— 8i?", we have, due to 



A, 



0: 



±2Vl^l - ^ dla r > C, 
" ' 2VI^1 + ^ dla r < C, 



I' 



(3.11) 



where the sign ± appearing in the flrst line depends upon 
the sign of Outside of the shell this sign may change at 
points where expression under the square root vanishes. 
On the contrary, inside the shell this sign is always posi- 
tive, because for a space-like hypersurface in Minkowski's 
space-time / is always increasing. 

The value of H is equal to the ADM energy calculated 
at r — > CX3. This is implied by the fact that, for large 
values of r, we have l{r) ~ r^. Consequently, due to 
p. lip , we obtain: 



1 



533 = n 



2 2H 



(3.12) 



We expect that the ADM mass H will play role of the 
Hamiltonian of the total "gravity + matter" system (cf. 
[3, |T^). Numerical value of H may be obtained from 
the singular part of equation (|2.16p . Since the singular 
part of the right hand side of this equation is equal to 
the jump — , we obtain: 



Bttc = 2R 




(3.13) 



where e denotes the sign of outside the shell. Now 
we insert the value of A_ calculated from (|3.5p (where 



5 



= 0), our "equation of state" (|2.12p and the value of 
p calculated from p.4p . This way we obtain: 



2H 



-C/ = e^l- — - v/l + C/2 , (3.14) 
and, consequently: 

[/) = I |l-(v/l + C/2 -[/)'! . (3.15) 
The value of e can be obtained from equation (|3.14p : 

e ^ sgn ( Vl + [/2 - [/ ) = 1 . (3.16) 

In the next Section we prove that the reduced phase 
space for gravitational degrees of freedom V may be 
parameterized globally by the two variables R and [/, 
whereas the function I and constants Cj P play role of 
gauge parameters. Given the point (i?, U) of this space, 
the specific choice of these parameters allows us to recon- 
struct the entire Cauchy data, with of A± and H given 
uniquely by equations (|3.6p and (|3.15p . However, gauge 
parameters are not completely arbitrary. One condition 
is obvious: Z(C) = B? . Moreover, inside the shell, I must 
increase monotonically from to B? . Because the sign 
"e" , given by p.l6p . is positive, outside of the shell I must 
increase and, due to asymptotical flatness, must behave 
like at infinity. Any function / is allowed provided it 
fulfills these conditions. Equations (|3.1ip . (|3.6p and p.ip 
enable us to reconstruct completely the data {n,l,f,u). 
Within our gauge subspace, given by condition (|3.ip . all 
states may be obtained this way. In the sequel we show 
that the states of the system obtained for other values of 
C and I, are equivalent up to a gauge transformation. 

For negative values of /3 equation p.6p implies that the 
entire external curvature vanishes outside of the shell: 
P*-'' = 0. Due to the fact that the external portion 
of space-time is Schwarzschild with mass iJ, the only 
surfaces satisfying this condition are standard surfaces 
{t = const}. Cauchy surfaces corresponding to different 
negative values of /3 coincide, therefore, outside of the 
shell and differ only inside the shell. 



IV. CANONICAL STRUCTURE OF THE 
REDUCED PHASE SPACE 

The dynamics of our system will be given uniquely by 
Hamiltonian p.isp . if only we know reduced symplec- 
tic structure expressed in terms of gravitational variables 
(i?, U) and, possibly, additional matter variables. 

It is possible that there is no matter at all, and matter 
Lagrangian is equal to zero. Yet, the space-time may 
still exhibit singularity and our "null shell" becomes a 
shock wave of pure gravitation. Phase space is then two- 
dimensional and symplectic form consists of gravitational 
part only. Its construction is given below. 



A. Reduction of gravitational part of the 
symplectic form 

We restrict gravitational part of the symplectic form 
il given by (|2.17p . to the gauge space (|3.ip and express 
it in terms of parameters {R,U,l,Q. It turns out that 
this form does not depend upon I and C- This proves 
that the latter are, indeed, gauge variables. For techni- 
cal reasons it is easier to work with the 1-form O given 
by equation (|2.18p . because the restriction to the gauge 
space commutates with exterior derivative of the form. 

Formula (|2.18p may be rewritten as follows: 



e 



n 



I 



AttR^PkSz 



K 



(4.1) 



Due to gauge condition (|3.ip . the first term in the above 
formula vanishes. From the equation (j3.6p we have: 



fl^ ^2UR^+^ B{C~r) , (4.2) 
where B denotes the Heaviside function. Then 

{S {fl^) } (r) = 2i?(C - r)S{UB'+^) , (4.3) 
and, consequently, 

(n/-^) {r)dr^6iUB^+^) , (4.4) 

because of the condition — B^ and due to vanishing 
of / outside of the shell. Finally we obtain 




e = -w6{UB^+^) + AttB'^PkSz 



K 



where w denotes 



{nl 2^ {r)dr 



Using equation (|3.1ip we obtain: 



- / 2 I'dr = 



/ 2 



(4.5) 



(4.6) 



I'dr 



(4.7) 

Because I is a monotonic function in interval [0, C] and 
— (H-/3) > 0, we may calculate the integral with respect 
to the following variable: 



'vr 

B 



-(1+/3) 



(4.8) 
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over the interval [0, ?7]. Hence 
1 



. (4.9) 



The last integral may be denoted as F{U), where F is 
the indefinite integral. It turns out that the specific form 
of F will not be needed. We have: 



1 



1 + 13 



F{U) {UR 



1+I3\i+B 



^F{U)5 [U^ 



(f({/)C/ttjj i?2^ - i R^U^^ F' {U)5{U) 

1-8 (f{U)U^ rA - I R^ , ^ <5(C/) 

(4.10) 



The first term in the above formula is a complete (vari- 
ational) derivative. Hence, it vanishes under the exte- 
rior derivative, when we calculate the symplectic form: 
— (58. The second term alone gives us, therefore, an 
equally good primitive form 8 for il: SQ = 58 = fl. 
Taking into account that 



1 



we obtain the following formula: 
1 



S{U) = 6 arsinh U 
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p Sn + AttpPkSz^ , 



(4.11) 



where by fi we denote the momentum canonically conju- 
gate to the variable p := R^: 



p := ar sinh U 
Consequently, we have: 

n = dQ=^6pA6p + Att5{pPk) A Sz^ 



(4.12) 



(4.13) 



Finally the Hamiltonian (|3.15p may also be expressed 
in terms of variables {p, p): 



H{p,p)^-^p{l-e-'^) 



(4.14) 



For the shock wave, when matter degrees of freedom 
vanish from the very beginning, the symplectic form is of 
the following form 



and Hamiltonian is given by formula (|4.14p . with p and 
p being now canonical variables. 



B. Geometric interpretation of the momentum p 

We are going to prove that the quantity p may be in- 
terpreted as an hyperbolic angle between the vector nor- 
mal to external Schwarzschild surfaces {tschw — const.} 
and the vector normal to internal Minkowski surfaces 
{^Mink = const.}. The angle q;(u, v) between two normal- 
ized vectors u,v is defined by their (hyperbolic) scalar 
product: 



coshQ;(u, v) :— (u|v) 



(4.16) 



Similarly as in Euclidean geometry, we call this quan- 
tity the angle between two surfaces: the Schwarzschild 
one and the Minkowski one. To prove this interpretation 
of p it is sufficient to use the formula ()3.1ip . On the 
internal side of the shell we use the second part of the 
formula and put I — R^ . We then obtain 



2^/I^i 



Vl + C/2 



(4.17) 



But inside of the shell the geometry of C is given by a 
three-dimensional spherically symmetric surface in the 
Minkowski space. It is a matter of a straightforward 
calculations to prove, that in Minkowski space-time the 
quantity on the left hand side of the above equation 
is equal to cosh a, where a is precisely the angle be- 
tween such a subspace and the Minkowski flat surface 
{^Mink = const.}. This implies that U = sinh a. But 
our surface is a smooth extension of the external side 
of Schwarzschild space {tschw = const.}. This finally 
proves that p = a \s the angle between leaves of three- 
dimensional foliations of Schwarzschild and Minkowski 
spaces. 



C. Reduction of material part of the symplectic 
form 



Let us first consider a simple non-trivial Lagrangian of 
the Dirac type in spinor theory (see Appendix [B|) : 



L = y^det gAsiz'^z^ — z^z^) . 
It implies second type constraints: 

Pi = z2 , P2 = -Z^ , 



that inserted into into symplectic form (I4.13P reduces it 
to the following form: 



n = ]^5p hSp + SirSi^z"^) A (5(7pz^) 



(4.15) Moreover, denoting 



P 

q 




(4.18) 



(4.19) 
(4.20) 
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leads to the following symplectic structure 



We also introduce new canonical variables 



ft — -Sfi A dp 



6p A Sq , 



(4.21) 



in a four-dimensional phase space parameterized by vari- 
ables: {p,q,^,p). Because Hamiltonian (|4.14p does not 
depend on the first pair of variables, hence their evolution 
is trivial: they remain constant in time. And evolution 
equations for "geometric variables" {fJ.,p) generated by 
Hamiltonian (|4.14p are of the following form: 



-2p, 



p = 2y/pe 



(4.22) 
(4.23) 



the same as in the case of a gravitational shock wave, 
described in the previous section. 

It turns out that a similar reduction is valid for a 
generic of Lagrangian dependent on two degrees of free- 
dom (Eq. (jA2[) ). but, typically, the momentum p canon- 
ically conjugated to p shall be modified by a function 
a{p, q), giving the new momentum u := p — a{p, q). Con- 
sequently, the Hamiltonian (|4.14p depends upon variables 
p and q via p = v + <j{p, q) and the dynamics of the ma- 
terial variables is no longer trivial. We can see that after 
performing reduction of material part of symplectic form 
(|4.13p for two degrees of freedom K = 1,2, it has the 
following form 



Vl = — A 5p 



Att5{pPk) a 5z 
From Eq. (jA2[) we have that: 

hence 



K 



(4.24) 



(4.25) 



1 



= inSp A Fa'(5z" -t- 'i'KpFK,LSz^ A + -Sp A 5p . 

(4.26) 

The assumption dF ^ implies that Fi \ — F1.2 7^ 
0. It is then either positive or negative, and the form 
Fk.lSz^ a 5z^ is a symplectic form in two dimensions. 
Besides, Darboux theorem implies that there exists such 
a coordinate system ^ ~ ^(z^,z^) and rj ~ rji^z^^z"^), in 
which 

F^f,i<5z^A(5z^ = SiFjiSz"^) = ^5 {^Stj - tj60 ■ (4.27) 

on 



Therefore 



FkSz"" = — {^Stj - r^60 + ^Sa , 
8n Stt 



(4.28) 



where C, = cr{£,,'ri) is a function of variables {S,,rj). In- 
serting the above two equations to the formula (|4.26p we 
thus obtain Vl of the following form: 



n = 6(^0 A <5(VP^/) + t:S{p -a) A dp 



(4.29) 



Vp V , 

p — a 



in order to simplify the form fi: 



Q = 5p A 5q+ -Su A Sp . 



(4.30) 
(4.31) 
(4.32) 



(4.33) 



Hamiltonian of the system, still equal to the same expres- 
sion (|4.14p . can be expressed in term of new variables as 
follows: 



H 



(.,p,p,g) = ivp(^l-e-^(''+K*'*))) . (4, 



34) 



This is universal form of Hamiltonian for spherically sym- 
metric self-gravitating shell of null matter coupled to two 
matter fields. Properties of any specific model of such a 
matter are uniquely implied by the function a, dependent 

upon two material variables This function is 

determined uniquely by the two functions Fk contained 
in the Lagrangian. 

Example: Consider the following Lagrangian density 
(see Appendix [B)) 



L = v/dct gAB {{z^ + 2zh^)z' + iiz^)^ ~ z')i^) . 

Rewriting this Lagrangian in new variables ^ and rj such 
that z^ = ^{S, + rj) and z^ = ^{^ — rj) we obtain 



The above equation expressed in terms of canonical vari- 
ables p = = ^{z^ + z2), q = = ^{z^ - z2) 
and p takes the following form: 

\o = \p''l\p + qW^q'). 
Hamiltonian of the whole system is following 



ii{v,p,p,q)^ ^■,Jp\\~ e 2p3/2 



V. DYNAMICS. RECONSTRUCTION OF 
SPACE-TIME GEOMETRY 



Hamiltonian (|4.34p uniquely generates dynamics of our 
system in the form of the following Hamilton equations 



8 



for canonical variables (p, and (p, q) : 



dp 



op ap 



(5.1) 

(5.2) 
(5.3) 



uq Oq 

The first two equations can be written in term of p = 
v+a and p only. Hence the equation (|5.1[) takes universal 
form, identical to (I4.22p : 

p = 2v/^e-2'^ . (5.5) 

Combining it with ([O]) and ((5^ leads to: 



da 



_9o_ 



da 
dp 



(5.6) 



Again, equations (|5.2p and (|5.6p imply the universal 
equation, identical to equation (|4.23p : 



/i 



1 1 



'1 



(5.7) 



We conclude that dynamics of "geometrical" variables 
{p, p) does not depend upon the choice of the model of 
matter and in the case p = v is identical to the one from 
Example. This dynamics is also identical to the case of 
gravitational shock wave in empty space-time, with no 
matter. 

Once we know the dynamics of variables (p, p) we can 
uniquely (up to the gauge) reconstruct the space-time, in 
which this dynamics is realized. Suppose that we have 
any explicit solution of equations (|5.5p and (|5.7p . Choose 
some gauge (3 < —1 and, separately for each moment of 
time, gauge variables Q and I. This enables us to recon- 
struct entirely the set of Cauchy data at each instant of 
time separately. To reconstruct the whole geometry of 
space-time we also need lapse and shift functions. For 
this purpose let us write Einstein equations in terms of 
canonical variables gki and P^^ . Because those objects 
are already known (as well as their time derivatives) in 
each moment of time, in result we obtain elliptic equa- 
tions for the lapse and the shift. In this way we obtain for 
the lapse function the second order equation in variable 
r as a condition for preserving of the /?-gauge in time. 
This equation has to be solved with the following bound- 
ary conditions: iV = 1 at infinity and ^ = at r = 0. 
In order to calculate the shift function we have to use the 
equation for time derivative of three-dimensional metric. 
This is the equation of the first order with respect to the 
shift function and enables us to reconstruct it uniquely. 



A. Solution of Hamilton equations 



From equations (|5.5p and (|5.7p : 



P = 2^e" 
1 1 



-2^^ 



{1 



there follows 



dp e^^M 
Solution of the above equation is of the form 



1 



1 - e^^Ai 



(5.8) 
(5.9) 

(5.10) 

(5.11) 



Therefore we have the equation for the time evolution of 
p: 



1 1 



2 

whose solution reads as 
1 



1 



log(e2^ - 1) - 



{t - to) 



(5.12) 



(5.13) 



The constant po may be expressed in terms of the whole 
energy of the system: y/po = 2E. Solution of equation 
(|5.8p for time evolution of p is following 



Elogp = t-to 



(5.14) 



TRANSITION FROM SCHWARZSCHILD 
TIME TO MINKOWSKI TIME 



VI. 



Equations (|5.5p and (|5.7p are the same as Hamilton 
equations for variables (/i, p) in case when a — 0. Then 
V — p and (/i,p) are canonical variables. Because of the 
identical form of dynamical equations, in further calcu- 
lations we limit to this case. 

Up to now, we have described the evolution of our 
system with respect to the Schwarzschild time, which 
coincides with the Minkowski time at spatial infinity. We 
can also describe our evolution with respect to any other 
time variable, corresponding to a different (3-|-l)-foliation 
of space-time. Leaves S'f of the new foliation does not 
have to coincide with previous hypersurfaces St — {t = 
const.}, but we assume that asymptotically (at space- 
infinity) they do. 

Change of the time variable is not a standard transfor- 
mation in classical mechanics. Here, we use an approach 
introduced by one of us in paper (llj , based on the notion 
of a contact manifold. 

Denote by v the retardation of the new time variable 
with respect to the previous one, calculated on the shell. 
This means that the hypersurface S^ intersects the shell 
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at Schwarzschild time t + v{p,ii). We assume that the 
value of this retardation depends upon the actual dy- 
namical situation, i.e. upon position and momentum, 
but it does not depend explicitly upon the time vari- 
able. This means that the new gauge conditions are 
intrinsic, depending only on initial data. The function 
V = v{p,p) contains the entire information about the 
transition between the old variables {p{t),fi{t)) and the 
new ones {p^ {t) , p^ (t)) , because once we know {p{t), p{t)) 
we can solve equations generating dynamics of a system 
and put 



p^{t) ■.=p{t + v{p{t),p{t))) 
pHt) ■.=p{t + v{p{t),p{t))) 



(6.1) 
(6.2) 



For an arbitrary, not constant, function w, such a trans- 
formation is, in general, not canonical. We will show in 
the sequel how to specify the canonical structure of our 
reduced phase space in terms of these new variables. For 
that purpose it will be convenient to use the language 
of contact manifolds. Observe that the entire informa- 
tion about dynamics of a system may be retrieved from 
the three-dimensional contact space, defined as a surface 
{E = H{p,p)} in the four-dimensional space {t^ E, p, p} 
equipped with the standard contact form: 



1 



-Sp h5p - 6E hSt 



(6.3) 



This symplectic form in four-dimensional phase-space be- 
comes degenerate when restricted to the surface {E ~ 
H{p, p)}. Trajectories of the system are defined uniquely 
as those whose tangent vector belongs to this degeneracy. 
To prove this it is sufficient to parameterize our subspace 
in terms of three variables {t, p, p), and rewrite the form 
in the following way: 

■.= -SpASp - l—Sp+—Sp] A5t . (6.4) 

We see that the vector annihilating the above form must 
be proportional to the vector 



^ d . d . d 
dt ^ dp ' dp 



(6.5) 



where p and p are given by Hamilton equations (j5.5p and 

(EZD. 

In a paper [Tl| it was shown that the choice of new 
gauge condition Q is equivalent to the choice of variable 
T 



T ■.= t-v 



(6.6) 



as new time variable. Let us rewrite, therefore, our sym- 
plectic form 4" in terms of the new time T. It will be 
convenient to have energy E instead of p as independent 
parameter and treat p as a, function p — p{E, p) obtained 



from solving equation (|4.14[) . Therefore we have that 



1 / dv dv 

* ^-5p ASp- SE A6T - SE A { —SE + —Sp 

2 \oE op 

=-(sp- 2^Sp] ASp-SE AST . 
2 V dp J 



(6.7) 



Now define: 



/dv 
— {E,p)dp + a{p) 



1 dv 



RdR 

where a{p) is arbitrary. Then, we have 
■ii^ypASp-SEA ST 



{E, R)dp + a{p) , (6.. 



where 



p := p -V {H{p,p),p) 



(6.9) 



(6.10) 



plays role of the momentum canonically conjugate to 
p^{t). Differentiating (|6.10p with respect to E provides 

4(A-.)^-2|(i^,p) = -^|^(i^,i^). (6.11) 

It is interesting to notice (cf. [ll|) that any choice of a 
function p — p{p, p) leads to a certain gauge condition, 
i.e. to the definition of a new time parameter on the 
shell, because we may always reconstruct the retardation 
function v from (|6.1ip . As an example consider p — 
Un{() = Airp/R (see formula (|3.4p ). which corresponds to 
the Minkowski time calculated at the internal side of the 
shell. Indeed, in Minkowski spacetime the momentum 
canonically conjugate to the shell's position R is equal 
to the kinetic momentum p. Change of variables from 
R to p implies transformation of momenta, according to 
identity: 



pdR 



2R 



dp , 



which, finally, gives ■^pdp when integrated over the shell. 
We are going to show in the sequel, that the evolution 
defined this way corresponds, indeed, to Minkowski time. 
For this purpose we express the new momentum in terms 
of the old one using (|4.12p : 



/ 2E 1 
p = sinhpJ 1 — = sinh/i e~^' ^ 2 ^"'^ ^ e"^'') , 

(6.12) 

where the energy E is equal to the value of the Hamilto- 
nian H. But, from (|4.14p . we obtain: 



2E 



(6.13) 
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and, consequently, 



E 



dji _ 1 
'dE~R 



On the other hand, we get from (|6.13p : 

dfi_ _ 1 
dE~ {R- 2E) ■ 

Thus, equation (|6.1ip imphes: 



2E 



dv 



dR R-2E 



.14) 



(6.15) 



(6.16) 



To prove that the above quantity describes retardation 
between the Minkowski time T inside the shell an the 
extrinsic Schwarzschild time T, differentiate the equation 
(|6.6p over t: 



dv ■ 

^-'-OR^ 



(6.17) 



where we have used the energy conservation E = 0. 
Derivative R can be calculated from equation (|5.5p . and 
relation l|6.13p . leading to: 



and finally 



dv ■ _2E 
OR IT 



2E 

T = 1 

R 



(6.18) 



(6.19) 



Integrating equation (|6.16p we get retardation between 
the Minkowski time T and the Schwarzschild time t: 

v{E, R) = 2E\og{R - 2E) , (6.20) 

hence 

T = t+2Elog{R-2E) = t+2RfilogR{l - 2fi) , (6.21) 

and R > 2E, otherwise /i would not be defined (Eq. 
I6.12p . Hamiltonian may by expressed in terms of 
Minkowski variables (p, /i): 



and Hamilton equations take the form: 

l_d_ 
2dT' 
1 d 



-P = \fp . 

1 fL 



2dT^ 2^ 



(6.22) 

(6.23) 
(6.24) 



Using relations = -^^i ^"^^ (|6.2ip between t and T 



we can reconstruct equations (15^5)1 and ([5? 

Dynamical equations (|6.23p and (|6.24p may be explic- 
itly solved: 



= r - To + 
P = Mo 



T-To + v^A) 



(6.25) 
(6.26) 



Keeping in mind that ,/p = R, the above solution for 
R (linear propagation with velocity equal 1) is implied 
by the fact that the only null-like, spherically symmetric 
surfaces in Minkowski space are light-cones. Hence, the 
evolution cannot be global: it ends at the cone's vertex. 
Interesting case of the global (in Minkowski time) evolu- 
tion is obtained when the two cones: the future oriented 
one and the past oriented one, cross each other. The the- 
ory of crossing shells was thoroughly analyzed in paper 



APPENDIX A: 



LAGRANGIAN FOR NULL 
MATTER 



Because the shell metric gab, a,b — 0,1,2; is de- 
generate, we have a priori no standard scalar density 
on S which can be used in the definition of the La- 
grangian. We must, therefore to manufacture such a 
density from the following ingredients: matter fields 
and their derivatives along shell, the degeneracy 
field X of gab on S and two-dimensional volume form 
A := \/detgAB on each surface = const.}. Choose 
the field X (otherwise given up to a multiplicative con- 
stant) in such a way that < dx^\X >= 1. Then, the 
following function: L = )^X°'z^, {K = l...n), 
is a scalar density and may be taken as a matter la- 
grangian. A more sophisticated example is given by 

L = x(^X''z^aX''z^f,aKL{z^)y {X'^^^cf, where ^ is a 

scalar field, 2a + P = 1 and (Tkl is a metric tensor in 
the space of field values. We may even generalize those 
examples to products of them. However in case of spher- 
ical symmetry we have X° = Sq and all of them give the 
trivial Hamiltonian H — 0, what implies equation A = 0. 
The shell surface is, therefore, an isolated horizon (cf. 
M, Eq. (6.2)) 

We are interested in examples of a null-matter, which 
may couple non-trivially to a generic null shell, not nec- 
essarily an isolated horizon. Such a nontrivial example 
can be obtained if one considers a thin shell of matter 
described by the Lagrangian density with at least two 
degrees of freedom, which depends not only on veloci- 
ties, but also on the configuration variables. Consider 
the Lagrangian density of the type: 



L = L{z 



\x 



aK 



(Al) 



where K, L = 2 . . . n, and Fk is a covector field defined on 
the space of material variables. We show in a sequel that 
this leads to a non-trivial model already for two degrees 
of freedom: K, L = 1,2. In spherically symmetric case 
this Lagrangian takes the form: 

L = \z'^Fk{z^) ■ (A2) 

Euler-Lagrange equations for this system read as 

AFi = A(F2.i - F^a)z^ , (A3) 

AF2 = -A(F2,i - Ti,2)i^ . (A4) 
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If the quantity (-F2_i — Fi^2) is equal to zero, i.e. if F is a 
closed form: dF = 0, then A = 0, and matter described 
by this Lagrangian density couples again to an isolated 
horizon. Assume that dF ^ 0. Then the Euler-Lagrange 
equations imply the following constraints equations: 

i^ii^ + Fi-P- = . 

□ 

APPENDIX B: EXAMPLES OF NON-TRIVIAL 
LAGRANGIANS 

1. A simple example of Lagrangian of the type (|A1[) 
may be the following one: 

L ^ \{z'^ - 7?) , (Bl) 

whose structure resembles the Dirac Lagrangian for 
spinor fields. This Lagrangian implies the following con- 
straints: 

AP2 = P2 = = -Azi (B3) 



These are second type constraints. Inserting them into 
the symplectic form ()4.13p we obtain: 



n=-S^iA5p + 87:6{y/pz^) A 5{^z^) . (B4) 



Similarly as in spinor theory, one configuration variable 
(in electrodynamics it is the imaginary part of spinor 
in Majorana representation) becomes momentum canon- 
ically conjugated to the second variable (the real part in 
Majorana representation). 

2. Another example: 



L = \{{z^ + 2z^z^)z^ + ((z^f - z^)i^) . (B5) 



is used to reduce the symplectic form described in Section 
HVCl 
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